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Based on density functional theory the noncentrosymmetric superconductor M03AI2C in its well 
established /J-Mn type (P4i32) crystal structure is investigated. In particular, its thermodynamical 
and dynamical stabilities are studied by calculating lattice vibrations within the harmonic approx- 
imation. It is found that the fully stoichiometric compound is dynamically unstable. However, 
compounds with carbon vacancies, i.e., MosA^Ci-^, can be dynamically stabilized for vacancy 
concentrations x > 0.09. By means of a simple thermodynamical model we estimate x ~ 0.13 — 0.14 
for Mo3Al2Ci_a; at the experimental preparation temperatures. The influence of the carbon vacancy 
concentration on the electronic structure is investigated. 
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I. INTRODUCTION 

M03AI2C has been already synthesized in the 19603^ 
and subsequently classified as a superconductor Re- 
cently it has attracted renewed attention because 
its cubic /3-Mn type crystal structure does not contain 
a center of inversion. In such noncentrosymmetric super- 
conductors, as first discussed for CePtsSip' the electron 
pairing process is described to be be a mixture of spin- 
singlet and spin-triplet statesP^E^ The question whether 
M03AI2C can be classified as a conventional or unconven- 
tional superconductor remains still unresolved according 
to very recent investigations!'^ : ! In this work, we will not 
tackle this issue directly but will provide results of ex- 
tensive density functional theory (DFT) calculations on 
the thermodynamical and dynamical stability as well as 
the electronic structure of this compound as a function 
of carbon content. 

M03AI2C crystalizes in the cubic /3-Mn type P4i32 
structure containing 24 atoms in the unit cell with, 
namely 12 Mo, 8 Al, and 4 C atoms. The C atoms are in 
the center of regular M06 octahedrons which are tilted to 
each other.™ The comparison of our DFT structural pa- 
rameters to experimental values^ in Table [I] shows both 
in excellent agreement. 

II. COMPUTATIONAL DETAILS 

The DFT calculations were done using the Vienna 
ab initio simulation package (VASPp 4 ^' utilizing the 
pseudopotential construction according to the projector 
augmented wave method.^ For the exchange-correlation 
functional the generalized gradient approximation as 
parametrized by Perdew, Burke, and Ernzerhofl^ was 
chosen. The potential for Mo contained 9 valence states 
including the 4p semicore states, whereas for Al and C 
three and four valence states were included, respectively 
The size of the basis set was defined by an energy cutoff 
of 500 eV. The Brillouin-zone integration for the com- 



TABLE I: Structural parameters and Wyckoff positions of 
M03AI2C. 
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"experimental results from Ref. \E\ 
'DFT calculation (present work) 



putation of total energies was made using the tetrahe- 
dron method with with Blochl's corrections^ based on a 
13 x 13 x 13 Monkhorst and PacH^ fc-point mesh, whereas 
for the structural relaxations and for the derivation of the 
force constants the first order Methfessel-Paxton smear- 
ing methocPSl on a 7 x 7 x 7 fc-point mesh was chosen. 

The vibrational properties were calculated within the 
harmonic approximation by making use of the direct 
force-constant method as implemented in our program 
package fPHON (full-symmetry PHON), which is based 
on the package PHONpHThe structural parameters, i.e., 
the volume and shape of the unit cell as well as the posi- 
tions of the atoms within the unit cell, were relaxed until 
the residual forces were less than 10 -4 eV/A. Further- 
more, for a high accuracy of the phonon calculations the 
force constants derived from the displaced atoms were 
corrected by subtracting the tiny but still finite forces 
of the fully relaxed structures. Some of the phonon dis- 
persions were cross checked by using density functional 
perturbation theory^ (DFPT) as implemented in VASP. 



2 



III. VIBRATIONAL PROPERTIES 

It turns out that perfectly stoichiometric M03AI2C is 
dynamically unstable. According to panel (a) of Fig. [I] 
optical modes with imaginary frequencies around T arise. 
The dynamical instability seems surprising considering 
the well-established crystal structure of Mo 3 Al 2 CPE2l 
One should, however, be aware that the actual carbon 
content is experimentally difficult to discern by X-ray 
diffraction because of carbons comparatively small X- 
ray cross section^ and therefore there exists some un- 
certainty with respect to carbon vacancies. In general, 
many carbides are prone to have vacancies on the C sub- 
lattice. 

Before speculating on the physical explanation of the 
detected instability we undertook numerical and method- 
ological tests. First of all, it is assured that the re- 
sults are converged with respect to the number of k- 
points. Furthermore, for deriving the force constants 
several atomic displacements of the atoms, i.e., |u| = 
{0.01,0.02,0.05,0.1,0.2,0.3,0.4} A were chosen and cal- 
culations without any symmetry were performed. All 
calculations confirm the existence of imaginary modes for 
M03AI2C. In particular, the calculation with the smallest 
displacement of |w| — 0.01 A yielded imaginary optical 
modes with a frequency of 1.62 i THz at T. As a further 
test the DFPT^22l technique as implemented in VASP was 
applied also resulting in imaginary modes of 1.65 i THz at 
r. From these DFPT calculations force constants have 
been derived and used as input for /PHON. Both, the 
direct force-constant method and the DFPT treatment 
gave very similar results for the phonon dispersions as 
shown in panels (a) of Fig. [I] 

In order to study if another similar structure exists 
that is dynamically stable and energetically more favor- 
able compared to the cubic /3-Mn type structure, all 
atoms were displaced from their equilibrium positions. 
In addition, a tetragonal deformation was enforced onto 
the unit cell. Using VASP this deformed structure was 
subsequently relaxed without any symmetry constrains, 
resulting in the well-known crystal structure described 
before. 

Because the perfectly stoichiometric compound 
M03AI2C is found to be dynamically unstable, we 
investigated if dynamical stabilization can be achieved 
by vacancies, in particular by carbon vacancies. 

In general, vacancies on all sublattices will stabilize 
the phonon dispersion, at least above a certain concen- 
tration. Vacancies on the Mo or Al sublattice are less 
likely to exist since they would be easily detectable, i.e., 
they have a comparatively large X-ray cross section!^ 
Furthermore, our DFT derived vacancy formation ener- 
gies in the following Section [TV] strongly indicate that Mo 
or Al vacancies are thermodynamically too costly to be 
formed. 

Assuming a certain carbon vacancy concentration suit- 
able supercell calculations were done for the defect struc- 
tures. Panel (b) in Fig. fTlshows the phonon dispersion for 



M03AI2C0.75, i.e., for a single carbon vacancy in the stan- 
dard unit cell (lxlxl supercell) with 3 out of possible 4 
carbon sites occupied. Further calculations, of which the 
dispersions are not shown, were also done for a single car- 
bon vacancy in a 2 x 1 x 1 supercell (M03AI2C0.875 with 
7 out of possible 8 C sites occupied) and in a 2 x 2 x 2 
supercell (M03AI2C0. 96875 with 31 out of 32 C sites occu- 
pied). Both M03AI2C0.75 and M03AI2C0.875 are found to 
be dynamically stable with no imaginary modes, whereas 
M03AI2 Co. 96875 is found to be significantly unstable with 
the value of the lowest imaginary optical modes at T be- 
ing 1.26 i THz quite close to that of the perfect stoichio- 
metric compound (1.62 i THz). 

In Fig. [2] the normalized phonon density of 
states (DOS) of the dynamically unstable compounds 
M03AI2C1 and M03AI2C0. 96875 are compared to those 
of the dynamically stable compounds M03AI2C0.875 and 
Mo3Al2Co.75- While at low frequencies no Debye-like u> 2 
behavior is observed for M03AI2C1 and M03AI2 Co. 96875 , 
it is seen for the other two cases. For M03AI2C0.875 the 
Debye-like behavior is observed only in a rather nar- 
row frequency range up to 0.5 THz due to the soften- 
ing of the acoustical and optical modes. However, for 
M03AI2C0.75 the Debye-feature reaches up to 1.4 THz. 
The partial DOS reveals that Mo, as the heaviest atomic 
species, dominates the lower frequency spectrum (up to 
7.5 THz), whereas C, being the lightest atom, has con- 
tributions only at frequencies above 13.5 THz. Further- 
more, the carbon dominated frequency modes are shifted 
down by the introduction of vacancies from above 15 THz 
for Mo 3 Al 2 Ci to about 13.5 THz for M03AI2C0.75. Strik- 
ingly, in all the DOS's a pronounced Al peak at rj 
12.2 THz occurs. The Al spectrum is rather broad with 
significant contributions from 6 THz to 12.2 THz. Even 
below this range a telling small contribution is found in- 
dicating a hybridization with Mo modes. 

The results indicate, that there is a critical concen- 
tration of carbon vacancies below which the compound 
becomes dynamically unstable. Assuming that the fre- 
quency of the lowest optical mode at T scales linearly 
with the carbon vacancy concentration x, the critical car- 
bon vacancy concentration is estimated to be x CT a ~ 0.09 
(using the values of MosA^Ci-^ at x = 0, 0.03125, and 
0.125 as input; see Fig. [3]). 

The cause of the stabilization of the optical low- 
frequency Mo modes by the carbon vacancies is the 
changed Mo-C bonding in the Mo 6 C subunits. We ex- 
amined the relaxations occurring in M03AI2C0.75, i.e., 
when a single C atom is removed from one of the four 
MoeC subunits in the unit cell of M03AI2C1. Thereby 
a strong influence on the Mo-C bonds in all three re- 
maining MoeC subunits is observed because they share a 
common Mo atom with the defect subunit. By relaxing 
the atomic positions the C atom drifts into an off-center 
position within the remaining MogC subunits, increasing 
the average Mo-C bond length by 1.4 % and the corre- 
sponding octahedral volume by 3.6 %. The distortion 
of the remaining MoeC subunits seems to be the stabi- 
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FIG. 1: (Color online) Phonon dispersions up to 5 THz of stoichiometric M03AI2C (a) and of M03AI2C0.75 (b). Imaginary 
frequencies are shown as negative values. For M03AI2C in panel (a) the calculated phonon dispersion derived from the direct 
force-constant method (black lines) is compared to results from DFPT theory (red lines). For M03AI2C0.75 in panel (b) 
the dispersion relation as calculated from the direct force-constants method (black lines) is presented. It should be noted 
that because of the reduced symmetry due to the carbon vacancy the dispersions along the paths A* = ±[£, — £, — f] and 
S* = {±[£, £, 0], ±[£, 0,£], ±[0, — £, £]} with < £ < n/a differ to the other symmetry related fc-directions refered by A and E, 
respectively. 



lizing factor for the vibrational modes. Concomitantly 
we notice this distortion for the Mo-Al bonding, e.g., the 
two distinctive nearest Mo-Al bonds with bond-lengths 
2.84 A and 2.95 A in the fully stoichiometric compound 
are distorted in M03AI2C0.75 to lengths in the range of 
2.79 — 3.01 A. Such a distortion for the Mo-Al bonds has 
recently been indicated experimentally.^ 



IV. CARBON VACANCIES 

As discussed, a certain amount of carbon vacancies is 
needed to stabilize the imaginary optical modes. The key 
question is if the formation of vacancies is at all thermo- 
dynamically possible. This question is investigated by 
calculating vacancy formation energies and by means of 
a model. 

Within a standard DFT approach, the vacancy forma- 
tion energy e^ ac per atom X is defined by subtracting 
the total energy E'dft(Moi 2 A1 8 C 4 ) of the stoichiometric 
compound from the total energy Edft^Mo^AIs^ — X) 
of the compound with a vacancy of atom type X and 
adding the ground-state energy Edft(X) of the removed 
atom X by 

e* c = £ DFT (Mo 12 Al 8 C 4 -A) + £ DFT (A)- 

-£ DFT (Moi 2 Al 8 C 4 ) . (1) 

Because standard DFT calculations are strictly valid only 
at T = K no temperature dependency has been yet in- 
troduced. This can be done by considering the temper- 
ature dependent vibrational free energies F v ^o an d defin- 
ing the vibrational vacancy formation energy per X atom 



TABLE II: Vacancy formation energies in eV for M03AI2C 
for lxlxl and 2x1x1 supercells without vibrational 
contributions (DFT, T = K) and with / X (™ C )(T) at T = 
1523 K and 1773 K. 
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lxlxl 


2x1x1 




e 1.74 


0.86 


0.60 


0.56 


T 


= 1523 K: e + f(T) 1.58 


0.50 


0.20 


-0.39 


T 


= 1773 K: e + f (T) 1.51 


0.39 


0.08 


-0.66 


C 


vacancy concentration x for MosA^Ci-^: 


0.25 


0.125 



similar to Equ. [IJ 24 

/4(T) = J F vib (Mo 12 Al 8 C 4 -A) + F vib (A)- 

-^vib(Mo 12 Al 8 C 4 ) . (2) 

Both, s* ac and /^ C (T) are formulated for a standard unit 
cell. To derive results for smaller vacancy concentrations 
larger supercells are needed and the stoichiometrics in 
Equs. [T] and [2] have to be scaled accordingly. The refer- 
ence energies Euft(X) and F v ib(A") were derived for the 
ground states of body-centered cubic Mo, face-centered 
cubic Al, and C in the graphene structure. 

The vacancy formation energies in Table [XT] and Fig. [4] 
at T = K are strongly positive for all types of vacan- 
cies, whereby the Mo vacancy with its formation energy 
of 1.74 eV is by far the most unfavorable one. Carbon 
vacancies are the most favorable ones with a formation 
energy of 0.60 eV for Mo3Al 2 Co.75. This value is reduced 
by only 0.04 eV for the smaller vacancy concentration of 

Mo 3 Al 2 C . 8 75- 

The experimental samples were synthesized at T = 
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FIG. 2: (Color online) Total and partial phonon DOS of Mo 3 Al 2 Ci (a), Mo 3 Al 2 Co.96875 (b), M03AI2C0.875 (c), and M03AI2C0.75 
(d) on the Total phonon DOS (black solid line), partial DOS of Mo (red solid line), partial DOS of Al (blue dashed line), and 
partial DOS of C (green dotted line). In the insets the total DOS at low frequencies is compared to a Debye-like uj 2 behavior 
(red dashed line). 



1773 K and heat treated at 1523 KP Therefore, theory 
needs to consider temperature dependent vacancy for- 
mation energies /^ C (T) combined with the composition 
dependent configurational entropy S'conf (x^P^ in order to 
compare with experiment. For the actual calculation of 
the vibrational free energy the small amount of imaginary 
modes in the fully stoichiometric compound M03AI2C1 
were omitted. 

Table [TT] and Fig. [4] show that the vibrational contri- 
butions reduce the strongly positive vacancy formation 
energies at T = K. While this reduction is rather small 
for the Mo vacancy (from 1.74 eV to 1.51 eV at 1773 K), 
it is much larger for the other two types of vacancies. 

In particular, the formation energy of the carbon 
vacancy in M03AI2C0.75 decreases from 0.60 eV to 
0.08 eV at 1773 K. Remarkably, this reduction is much 



larger for the smaller carbon vacancy concentration, i.e., 
M03AI2C0.875J with a decrease by more than 1 eV down 
to —0.66 eV. These negative values for the formation en- 
ergy indicate a possible thermodynamical stabilization of 
carbon vacancies in MosAlaCi-a;. Noticeably, this differ- 
ence of the temperature dependent free energies for dif- 
ferent carbon vacancy concentrations comes exclusively 
from the vibrational contributions /(x, T), because -as 
mentioned above- at T — K the vacancy formation 
energies are almost equal. 

From an isolated defect modeP^ the temperature de- 
pendent equilibrium vacancy concentration x can be 
calculated. However, the vacancy formation energy 
e{x) + f(x,T) in the description of the internal energy, 
U(x,T) = (e(x) + f(x, T))x^ is strongly dependent on 
x. Hence, the isolated defect model cannot be applied 
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FIG. 3: (Color online) Frequency of the lowest optical mode 
at T versus carbon vacancy concentration x as calculated 
(red circles) and linearly interpolated (solid black line). The 
predicted critical carbon vacancy concentration is indicated 
(green dashed-dotted line). 
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FIG. 4: (Color online) Vacancy formation energy e + f(T) as 
the sum of the DFT e and vibrational f(T) formation energy 
versus temperature for a Mo vacancy (red solid line), a Al 
vacancy (blue dashed line), a C vacancy (green dotted line) 
in a unit cell (lxlxl supercell), and for a C vacancy in a 
2x1x1 supercell (green dashed-dotted line). 



directly, as the internal energy U{x,T) is not a linear 
function of x. In our case it is described as a quadratic 
function of x, U(x, T) — ax 2 + bx wherein a and b are 
temperature-dependent parameters fitted to the calcu- 
lated carbon values at x = 0, 0.125, and 0.25. 

Thus, the free energy for the vacancy formation includ- 
ing the configurational entropy is formulated as 

F(x, T) = ax 2 + bx — kBTS con f(x) . (3) 

Assuming S con f(x) is the mixing entropy of non- 
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FIG. 5: (Color online) Temperature dependent C vacancy 
concentration x(T) plotted as a solid red line. Experimen- 
tal preparation temperatures are indicated as dashed black 
lines. The critical C vacancy concentration, below which 
Mo3Al2Ci-a; gets dynamically unstable, is shown as a green 
dash-dotted line, while the C vacancy concentrations of the 
calculated supercells are drawn as blue dotted lines. 



interacting vacancies S con f(x) = a;ln(a;) + (l— x) ln(l— x), 
the derivative of F with respect to x can be used to 
search for the temperature dependent concentration x(T) 
by minimizing the free energy, i.e., 

dF(x,T) _ e-zwt) 

dx l + e -l3(2ax+b) ' W 

with P = l/(kBT). This expression enables the numeri- 
cal calculation of x(T), as shown in Fig. [5] using a brack- 
eting root finding algorithm.^ Inspecting Fig. [5] the C 
vacancy concentration is 0.13—0.14. at the experimen- 
tal preparation temperatures^ of 1773 K and 1523 K. As 
elaborated in the previous Section |III| at such vacancy 
concentrations MoaA^Ci-a, is dynamically stable. 



V. ELECTRONIC STRUCTURE 

After finding that vacancies do exist on the carbon sub- 
lattice and that these are necessary to stabilize the crystal 
structure we will now briefly discuss the band structure 
and the electronic DOS of M03 A^Ci-^ . As will be shown 
the influence of a changed carbon stoichiometry on the 
band structure can not be described by a simple rigid 
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FIG. 6: (Color online) Electronic band structure and DOS of stoichiometric M03AI2C (a) and of M03AI2C0.75 (b) calculated 
scalar relativistically (black lines) and fully relativistically including spin-orbit coupling (red lines). The Fermi energy _Bf 
referring to the number of valence electrons of M03AI2C is indicated as a dotted line while the Fermi energy corresponding to 
M03AI2C0.75 is indicated as a dashed line. Note the different directions E, E* and A, A* as discussed in the caption of Fig. [I] 



band model. Especially, the spin-orbit splitting on the 
bands in a fully relativistic calculation strongly depends 
on x. 

In Fig. [6] the electronic band structure and DOS of 
the fully stoichiometric compound is compared to that of 
Mo3Al 2 Co.75- The attentive reader might question this 
choice of x = 0.25 as being too high, because our ther- 
modynamical model (described in the previous section) 
predicted a much lower carbon vacancy concentration. 
However, we chose this value due to the fact that the 
unit cells of both compounds have equal size and shape 
making it easier to compare these. 

At this point it should be remarked that our calculated 
band structure for M03AI2C1 shown in panel (a) in Fig. 
[6] resembles the one published previously by Bauer et 
aiP but differs distinctively from that published by Karki 
et al^ The different finding by Karki et al. can only 
stem from the use of a different crystal structure. We 
have recalculated the band structure by means of the full- 
potential linearized a ugmen ted plane- wave method using 
our own code FLAIRpE! an d also by comparing to a 
recent calculation with the code Wien2KP^ (the later was 
also used by Karki et alty . All these calculations yielded 
the same result for the band structure of Mo 3 Al 2 Ci in 
the cubic /3-Mn type crystal structure, i.e., the one shown 
in panel (a) of Fig. [6] 

We have calculated the band structures both in a scalar 
relativistic approximation, omitting spin-orbit coupling, 
and fully relativistically, including spin-orbit coupling in 
a self-consistent manner. As expected, some degenera- 
cies at the high-symmetry points are different whether 
spin-orbit coupling is included or not, e.g., a splitting 
of 30 meV at V and of 53 meV at R in M03AI2C1 and 
of 20 meV at both T and R in M03AI2C0.75 just be- 
low the Fermi energy E-p is evident. However, the most 
striking point is the loss of the double degeneracy of the 
bands due to spin-o rbit coupling in noncentrosymmet- 
ric compounds ! 31 * 32 ! From Fig. M we observe this vertical 



spin-orbit splitting of the bands, e.g., 65 meV on the path 
A around —0.35 eV in Mo 3 Al 2 Ci and 90 meV on the path 
A* around —0.3 eV in Mo3Al 2 Co.75- As a consequence 
the Fermi surfaces of noncentrosymmetric Mo3Al 2 Ci_a; 
do also split due to spin-orbit coupling which can be seen 
in Fig. [6] as the horizontal band splitting at Ep ■ 

Comparing the band structure and electronic DOS of 
Mo3Al 2 Ci with that of Mo3Al 2 Co.75 in Fig. [6] one im- 
mediately notices that removing a carbon atom from the 
carbon sublattice has a substantial effect and it is not suf- 
ficient to simply shift the Fermi level in order to account 
for the carbon vacancy. 

From these results we conclude, that the structure of 
the Fermi surfaces and hence any nesting, paramount to 
the understanding of superconductivity, strongly depend 
on x. 



VI. CONCLUSIONS 

In the present work we have shown, that vacancies are 
necessary to dynamically stabilize the cubic /3-Mn type 
(P4i32) crystal structure of Mo 3 Al 2 C, whereby vacancies 
on the carbon sublattice are energetically the most fa- 
vorable ones. According to our thermodynamical model 
the most probable carbon vacancy concentration x in 
MogA^Ci-a is about 0.13 — 0.14 considering actual ex- 
perimental preparation temperatures. 

We have demonstrated that there exists a critical value 
of x cr i t ~ 0.09 below which Mo3Al 2 Ci_ x becomes dy- 
namically unstable, and especially the frequency at which 
the Debye-like behavior of the phonons ends strongly de- 
pends on x. 

Likewise, the band structure and electronic DOS are 
influenced by the carbon vacancy concentration. 

The still unresolved question if M03AL2C is a conven- 
tional or unconventional superconductor can only be an- 
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swered when the carbon vacancies are properly consid- 
ered, as the structure and nesting of the Fermi surfaces 
depend on the carbon vacancy concentration. If this car- 
bon vacancy concentration could be controlled, it might 
be possible to tune the superconducting properties of 
MoaAlaCi-a;. This might be rather difficult, because 
it may depend on the sample preparation and the cool- 
ing process. Further, one can safely assume, that at the 
preparation temperature of ss 1500 K the sample is in its 
thermodynamic equilibrium. This is not the case when 



its superconducting properties are measured where it is 
in a quenched meta-stable state. 
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